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ABSTRACT 

We simulate the magnetosphere of the nearby millisecond pulsar PSR J0437- 
4715, which is expected to have an unscreened electric potential due to the lack of 
magnetic pair production. We incorporate General Relativistic (GR) effects and 
study curvature radiation (CR) by primary electrons, but neglect inverse Comp- 
ton (IC) scattering of thermal X-ray photons by these electrons. We find that the 
CR spectrum cuts off at energies below ~ 17 GeV, well below the threshold of the 
H.E.S.S. telescope (< 100 GeV), while other models predict a much higher cutoff 
of > 100 GeV. GR theory also predicts a relatively narrow pulse {(3° ~ 0.2 phase 
width) centered on the magnetic axis. EGRET observations above 100 MeV sig- 
nificantly constrain the application of the Muslimov & Harding (1997) model for 
7-ray production as a result of GR frame dragging, and ultimately its polar cap 
(PC) current and accelerating potential. Whereas the standard prediction of this 
pulsar's 7-ray luminosity due to GR frame dragging is ~ 10% of the spindown 
power, a non-detection by forthcoming H.E.S.S. observations will constrain it to 
< 0.3%, enforcing an even more severe revision of the accelerating electric field 
and PC current. 

Subject headings: stars: neutron — pulsars: individual(PSR J0437-4715) 



1. INTRODUCTION 

Several authors have included General Relativistic (GR) frame dragging in models of 
pulsar magnetospheric structure and associated radiation and transport processes, recognis- 
ing it to be a first order effect (see e.g. Muslimov & Harding 1997 (MH97); Dyks, Rudak, & 
Bulik 2001). 
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Usov (1983) was the first to suggest that the low magnetic field strengths of millisecond 
pulsars (MSPs) allow 7-rays up to at least 100 GeV to escape pair production. Most MSPs 
have (largely) unscreened electric fields due to the low optical depths of primary curvature 
7-rays for pair production in such low-B pulsar magnetospheres (Harding, Muslimov, & 
Zhang 2002, HMZ02). Radiation reaction limited curvature 7-rays up to about 100 GeV 
from MSPs have been predicted (HMZ02; Bulik, Rudak, k Dyks 2000, BRD00), making 
nearby MSPs such as PSR J0437-4715 (Johnston et al. 1993) attractive targets for ground- 
based 7-ray groups (BRD00; Venter 2004). The unscreened case offers a test for fundamental 
GR electrodynamical derivations of the polar cap (PC) current and potential, without having 
to invoke additional modifications such as pair formation fronts (Harding & Muslimov 1998, 
HM98) with associated slot gaps (Muslimov & Harding 2003) to explain most observations 
of canonical (high-B) pulsars. 

The use of an unscreened GR electric field (see section 2) for PSR J0437-4715 (implied by 
its relatively low spindown power - HMZ02) was justified a posteriori (see section 3). Several 
important parameters, most notably its mass and distance (Van Straten et al. 2001), are 
accurately known, making PSR J0437-4715 one of the closest pulsars to earth and probably 
much brighter and easier observable than other MSPs. Also, observations show that the 
radio and X-ray beams virtually coincide (Zavlin et al. 2002), implying that the observer 
sweeps through the approximate center of the PC (Manchester & Johnston 1995; Gil & 
Krawczyk 1997). 

In this paper, we investigate the effect of GR constraints on MSP spectral cutoffs, pulse 
profiles, integral flux and conversion efficiency of spindown power to 7-ray luminosity by 
simulating (using a finite element approach) radiative and transport processes which occur 
in a pulsar magnetosphere. 



2. THE UNSCREENED ELECTRIC FIELD AND RADIATIVE LOSSES 

We use the GR-corrected expressions for a static dipolar magnetic field (e.g. Muslimov 
& Tsygan 1992 (MT92); MH97) and curvature radius p c (e.g. HM98) of an oblique pulsar 
with magnetic moment /i = B Q R 3 /2 inclined at an angle x relative to the spin axis. The 
value of the surface magnetic field (at the pole), B , was solved for using (MH97) 



with E rot the spindown power, Q the angular speed, Cl the time-derivative thereof, I the 
moment of inertia, R the stellar radius, c the speed of light in vacuum and f(rj) defined by 
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eq. (25) of MT92. 

The effect of GR frame dragging on the charge density, electric potential and hence 
the magnitude of E\\ (the electric field component parallel to the local magnetic field lines) 
was carefully modelled for the unscreened case, since the optical depth for magnetic pair 
production above the PC is insignificantly small (see section 3). The 'near' and 'distant' 
cases for E\\ (when 77 ~ 1 and r) >> 1, with rj — r/R), coincide at different points for different 
pulsar parameters. We use the same framework as Harding, Muslimov and Tsygan (MT92; 
MH97; HM98), with all the symbols corresponding to their formalism. For the 'near' case, 

E \T = { 12«©oSi c o s X + Gs 2 Q 3 H (1)5(1) sin x cos 0} , (2) 
R 

with the vacuum potential $o = B QR 2 /c, compactness parameter k = el /MR 2 , e = 
2GM/Rc 2 , pulsar mass M, polar angle of last closed magnetic field line 6(77) = [(QRrj/cf (77)] x l 2 
and 6 = 6(1) (HM98). Furthermore, in eq. (2), 

= f^rj.^dlHJi (4) 
kfJ 2 (ki) L J 

with fcj and fcj the positive roots of the Bessel functions J and J 1; ^ = 9/Q(i]) the normalised 
polar angle, the magnetic azimuthal angle, and H (1)5(1) ~ 1. (For 7^ and 7,, see eq. [22] 
in HM98 and definitions following eq. [43] in MT92). Note that EZ eax (r) = 1) = as required 
by the boundary conditions and that E\\ scales linearly with radial distance 77 close to the 
stellar surface (derived from a Taylor expansion of eq. [3] and eq. [4] at 77 ~ 1). For the 'far' 
case (77 > Rpc/R), we use (HM98) 

E t - i 1 -£ 2 ) o {^cos X + ^e(77)//(77)5(77)esin X cos0| , (5) 
and for the corotating charge density p e , we use eq. (32) in MT92. 

The change in the energy of a primary electron is given (when only the dominating cur- 
vature radiation (CR) component is considered, neglecting inverse Compton (IC) scattering 
and synchrotron radiation) by 

dE n „ 2 (e 2 c\ . 

— = eP r cE\\ — (6) 



dt " 3\p 2 c 

with e the electron charge, (3 r = v e /c ~ 1 the normalised electron speed and 7 the Lorentz 
factor. The photon energy e 7 is set equal to the characteristic CR energy e c = 1.5(A c /p c )7 3 
(in units of m e c 2 - Luo, Shibata, & Melrose (2000)), with A c = h/m e c ps 3.86 x 10 -11 cm the 
Compton wavelength. 
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3. PAIR PRODUCTION, SPECTRA AND CUTOFFS 

According to HMZ02, the CR death line is at E rot < 10 35 ergs -1 . Evaluating E TOt = 
—AnlP/P 3 ~ 4 x 10 33 ergs -1 (using the corrected intrinsic period derivative P - Van Straten 
et al. (2001)) suggests that no pair production will take place. Detailed modelling yields 
negligible optical depths, confirming this scenario. This is indeed fortunate because of the 
limited number of free parameters in this case. However, a low intensity of IC scattered 
UV photons/soft X-rays into the TeV range may contribute to a weak pair production 
component. 

For the parameter ranges R % = i?/10 6 cm = 1.3 — 1.7 (e.g. Kargaltsev, Pavlov, & 
Romani 2004), 7 45 = //10 45 g.cm 2 = 1- 3 (e.g. HMZ02), and ( x ,() = (35°,40°) (Manchester 
& Johnston 1995), ( X ,C) = (20°,25°) (e.g. Pavlov & Zavlin 1997) and ( x ,0 = (20°, 16°) (Gil 
& Krawczyk 1997), with ( the observer angle, the maximum CR cutoff energy is obtained 
by using i? 6 = 1.3, J45 = 3 and \ = 20°. We used M = 1.58M derived from Shapiro 
delays (Van Straten et al. 2001). This corresponds to k ~ 0.2 and surface magnetic field 
strength B 8 = B /10 8 G ~ 7.2 (see eq. [1]). The relat ive altitude for maximum CR energy 
is obtained as rj ~ 1.47 corresponding to a normalized field line colatitude of £ ~ 0.1 and 
p c ~ 10 8 cm, while the magnetic azimuth = results in a maximum GR potential. The 
analytical expression for the maximum 7-ray energy is obtained by combining eq. (5), (6), 
and the expression for e 7 , giving 

W = (£) ^ A c (^) ^ P 1 / 2 < 17 GeV. (7) 

One of the most interesting predictions from MH97 is that the primary electron lumi- 
nosity is given by (assuming x ~ 0): 

£p X r5max ~ I^ 1 ~ K )^rot- (8) 

It is important to note that the electric potential and charge density were derived assuming 
that electrons leave the PC with a speed equal to c. Even if the stellar injection speed 
(3rc c, it can be shown that the electrons will become relativistic very close to the neutron 
star surface, making maximum electron energies virtually independent of the injection speed 
(also: A.K. Harding 2004, personal communication). The bolometric particle luminosity of 
a single PC will therefore be given by (MH97) 

L prhn = ac J \p e \®dS, (9) 

with $ the electric potential and dS the element of spherical surface cut by the last open 
field lines at radial distance r. Integrating over £ and <p, and letting rj — > oo, we obtain 
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(Venter 2004) 

rlx=o f 2 /3e #(l)[7r/2-e #(l)]\ . 2 1 

Vinvnax = |c 0S 2 X + ^ ^ _ K ) j 8111 * j ' ( 10 ) 

providing we adopt a value of = n/ 2 for distances rj > c/QR. This result reduces 

back to eq. (8) when x is se t equal to zero. We calculated the maximum efficiency of 
conversion of pulsar spindown power into particle luminosity L pr im,max for x = 20° and 
X = 35°, and obtained ~ 2 - 11% for PSR J0437-4715, for each PC (depending on R and I, 
using M = 1.58M ). We also obtained the bolometric photon luminosity L 7 using a finite 
element (particle tracing) approach and integrating numerically over all photon energies and 
field lines from the surface to the light cylinder: 

n27r/ pr=c/Cl \ 

\N(<f> R ,0 R ) / P^,9,r)dt)d^d9. (11) 
y Jr=R J 

Here P 1 is the CR photon power integrated over frequency and N = p e cdS/e the number of 
particles ejected per second from a PC surface patch dS at r = R centered at (<fi R , 9 R ). Since 
we cannot start with (5 R — 1 (i.e. infinite Lorentz factor), we assumed values close to 1 and 
found convergent photon luminosities of 2 — 9% of the spindown power (depending on R, I, X 
and 0, i.e. £ 7 /£ pr im,max ~ 1- This means that almost all particle luminosity is converted to 
photon luminosity as expected for strong radiation reaction. Radiation reaction, combined 
with further (weak) acceleration towards the light cylinder, result in a total residual electron 
power of ~ 1 — 2.5% of the spindown power at the light cylinder. 

It should be noted that the fundamental unscreened expression for E» (eq. [5]) changes 
sign along ~ 40% of the magnetic field lines originating at the PC. This field reversal is 
most dominant when <fi ~ ir, whereas no field field reversals occur for <fi ~ 0. Trapping of 
electrons may ensue at magnetic field lines along which the electric field reverses. We expect 
the system to reach a steady state as a result of the redistribution of charges along these 
field lines. These lines may become equipotential lines, or a reduced current may develop, 
resulting in the suppression of particle acceleration along them. This justifies our neglect 
of these field lines when calculating the pulse profiles, bolometric photon luminosity and 
integral flux. 

Figure 1 shows the pulse profiles for different observer angles (, for x — 35°. Maximum 
observed photon flux is obtained for ( ~ x an d for large values of cos0 (as in eq. [5]). The 
'dip' in light curves with ( > x near phase 0l/27t ~ 0.5 (where ~ n) might be due to the 
sign reversal of the electric field, because the magnetic field lines where this sign reversal 
occur, were ignored as noted above. 
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The differential photon power gLL 7 (0£,£, E)/d(pLd(dE per phase bin d(f>L, per observer 
angle bin d(, per energy bin dE, is obtained by inserting the product of the ratios of indicator 
functions I{4>Li4>l + d<f>L), HCX + ^C) an d I(E,E + dE) and their respective bin widths 
d<pL, d( and dE in the integrand of eq. (11). This allows us to compare the expected integral 
photon flux with EGRET upper limits above 100 MeV and 1 GeV (Fierro et al. 1995), as 
well as with forthcoming H.E.S.S. observations of this pulsar (Venter 2004). Note that the 
imaging threshold energy of H.E.S.S. is ~ 100 GeV (Hofmann 2001), although a non-imaging 
"pulsar trigger" for timing studies down to > 50 GeV can be employed for pulsar studies 
with H.E.S.S. (de Jager et al. 2001). 

The phase-averaged photon flux (as would be seen on a DC skymap) for a single PC 
may be calculated by 



F°J>E) = 



(3° 



d 2 An 



E 



2tt 



1 

~E' 



dL(<f>' L ,C,E') 
d<p' L dC'dE' 



d(f)' L dC'dE', 



(12) 



with distance d = 139 pc, (3° = A(j) L /2ir, A(f) L the pulse width in radians, Ail (> E) = 
sin(d(A(f) L the beaming solid angle, and d( taken arbitrarily small. Only one PC is expected 
to be seen, given the relative orientations of the magnetic axis and observer line-of-sight to 
the spin axis. The superscript 'o' will be used to indicate quantities applicable to an observer 

withCe (CC + O- 

The energy spectrum dL/dE due to CR is quite hard, resulting in a constant time- 
averaged integrated photon flux F^(>E), seen by the observer, as shown in figure 2 (e.g. 
curves (a) and (b)). The 100 MeV and 1 GeV EGRET flux upper limits from Fierro et al. 
(1995) are indicated by the squares on figure 2, which clearly constrain the flux band defined 
by (a) and (b). If we define an a priori phase interval of f3° ~ 0.2, centered on the radio pulse, 
and recalculate the EGRET flux upper limits from the factor five (= 1/(3°) reduced skymap 
background, we should get the even more constraining upper limits given by the diamonds 
in figure 2. We therefore have to revise the predicted fluxes for PSR J0437-4715 and we do 
so based on the following scaling argument: If we assume that the particle and hence 7-ray 
luminosity only scales with the spindown power and neutron star compactness, as in eq. (8) 
and eq. (10), i.e. the product of the current and voltage for such a pair-starved pulsar is a 
constant as predicted by eq. (8), we may scale the set of curves (a) through (c) (according 
to this condition of a constant photon luminosity L°) in terms of the limiting voltage and 
hence the cutoff energy to give F° 1 (>i?i) x -E cuto fr,i ~ E° 2 (>E 2 ) x -E C utofr,2 (for constant 
(3° and An , and energies E\ < E 2 ; E cntoS) i < E cntoS)2 )- In particular, when curve (c) is 
scaled according to £ C utoff,2 = A£ cuto fr,i, implying F° l2 (>E 2 ) ~ F^ ^E^/X, with A = 400, 
curve (d) is obtained, which no longer violates the revised EGRET upper limit at 1 GeV, but 
the cutoff energy then shifts up to ~ 1 TeV. Furthermore, if curve (d) is now translated so 
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that the energy cutoff also falls below the H.E.S.S. sensitivity curves, curve (e) is obtained, 
which would be consistent with both EGRET and H.E.S.S. (if the latter instrument does not 
detect this pulsar). Also shown is the flux band calculated for PSR J0437-4715 by BRDOO. 
Again, for a power-law photon spectrum with exponential cutoff, it can be shown that 
F°(>Ei)E cuto g ~ j3°L°/d 2 AQ° (assuming E\ <C ^cutoff, and F°(>E) having a flat slope due 
to CR). In order to constrain PSR J0437-4715's bolometric photon luminosity by forthcoming 
H.E.S.S. observations, we postulate that L 7 = aL°, where a = a(xX) 3> 1 is a geometrical 
factor correcting from the incremental luminosity corresponding to the observer's line-of- 
sight, to the total 7-ray luminosity of the pulsar. It then follows that L 7 ~ xF°(>Ei)E cutoS , 
with x(x, C) — o:d 2 2n sin (d(, which was found to be more or less constant for the same x 
and (. A non-detection by H.E.S.S., as implied by curve (e), leads to a 7-ray luminosity of 
< 0.003-E r ot- This value should be compared with the prediction of L 1 ~ 3 x 10~ 5 -E ro t given 
by Rudak & Dyks (1999) for a canonical pulsar with P = 1 ms and Bq = 10 9 G and with 
L 7 ~ 0.04i? rot predicted for pair-starved pulsars with off-beam geometry (using P « 5.76 ms 
and E rot ~ 4 x 10 33 erg.s" 1 - Muslimov & Harding (2004)). 

4. CONCLUSIONS 

CR cutoff energies for MSPs such as PSR J0437-4715 were predicted to be in the range 
50 — 100 GeV by HMZ02 and BRDOO, making proposals for ground-based telescopes with 
imaging thresholds near 100 GeV (e.g. H.E.S.S. (Hofmann 2001) and CANGAROO (Yoshida, 
Yoshikoshi, & Yuki 2002)) attractive. From the present GR theory, one would conclude 
that these telescopes may not be able to see the spectral tail corresponding to the intense 
primary CR component, since the hard primary CR spectrum does not extend to energies 
above ~ 20 GeV, as verified by both analytical and numerical (finite element) approaches. 
An IC component resulting from TeV electrons scattering the UV/soft X-rays from the 
surface of PSR J0437-4715 may however still be detectable, although this prediction by 
BRDOO should also be re-evaluated within a GR electrodynamical framework. However, it 
is quite obvious that the predicted time-averaged observer flux violates the EGRET upper 
limit at 100 MeV, impyling a revision of the existing theory. Forthcoming H.E.S.S. and 
future GLAST observations will help to constrain the 7-ray luminosity, and therefore the 
accelerating electric field. 

The authors would like to acknowledge useful discussions with A.K. Harding, B. Rudak 
and A. Konopelko. This publication is based upon work supported by the South African 
National Research Foundation under Grant number 2053475. 
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Fig. 1. — Photon luminosity (in relative units) vs. observer pulse phase (with phase 0.5 
corresponding either to <fi = or <fi = tt, depending on () for PSR J0437-4715 for different 
C (see legend). The following parameters were assumed (see text for references): P ~ 
5.76 ms (period), Rq = 1.3, J45 = 1, M = 1.58 M & and \ = 35°. The radio pulse at 
4.6 GHz (thick solid line - Manchester & Johnston (1995)) is superimposed for reference 
(see www.atnf.csiro.au/research/pulsar/psrcat). The "valleys" at observer phase ~ 0.5 of 
the lightcurves with ( > x are probably due to electric field sign reversals (FSR), since the 
magnetic field lines where these reversals occur, were ignored (see text for details). 




Fig. 2. — Observer time-averaged integral flux vs. threshold energy. Curve (a), for which 
R G = 1.7, 7 45 = l,x = 35 ° and C = 40° and curve (b) for which R 6 = 1.3, I 45 = 3, x = 20° 
and ( = 16°, define a "confidence band" wherein the integral flux is expected to lie according 
to the GR model discussed in this paper. Curve (c), for which i? 6 = 1.5, J45 = 2, x — 20°, 
and C = 16°, represents an intermediate curve. Curve (d) is curve (c) scaled with scale 
factor A = 400, while curve (e) is curve (d) shifted to the left (see text for details). The 
band with dot-dashed curves is that of BRD00 for PSR J0437-4715 for their 'Model A'. The 
squares represent EGRET integral flux upper limits (Fierro et al. 1995), while the diamonds 
represent these upper limits reduced by a factor v^5, appropriate for a beam with main pulse 
width of ~ 0.2. Also indicated are the H.E.S.S. sensitivities for 50 hours (Hinton 2004) and 
8 hours observation time, and the energy above which pair production is expected to take 
place (BRD00). 



